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Partial CMB maps: bias removal and optimal binning of 
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ABSTRACT 

We present a semi-analytical method to investigate the systematic effects and statisti- 
cal uncertainties of the calculated angular power spectrum when incomplete spherical 
maps are used. The computed power spectrum suffers in particular a loss of angular 
frequency resolution, which can be written as 5i T^/"fmaxi where jmax is the effective 
maximum extent of the partial spherical maps. We propose a correction algorithm to 
reduce systematic effects on the estimated Ce, as obtained from the partial map pro- 
jection on the spherical harmonic Y^^ basis. We have derived near optimal bands and 
weighting functions in ^-space for power spectrum calculation using small maps, and 
a correction algorithm for partially masked spherical maps that contain information 
on the angular correlations on all scales. 

Key words: methods: data analysis -methods:statistical -techniques:image process- 
ing cosmology: cosmic microwave background 



1 INTRODUCTION 

The measurement of the temperature and polarisation 
anisotropies of the Cosmic Microwave Background (CMB) 
radiation provides essential information for testing the cos- 
. mological models and determining their parameters. Most 
' of the statistical information present in the CMB tempera- 
, ture or polarisation sky maps can be encoded in the angular 
I power spectrum Ce- An overview of the physical mechanisms 
responsible for the CMB anisotropies and the effect of the 
cosmological parameters on the pow er spec trum shape can 
be fo und in (|Zaldarriaga et a]||l997l ') and (|Hu fc DodelsonI 
I2OO2D . 

The following references give an overview of some re- 
cent CMB power spectrum measurements: The all sky 
WMAP (Wilkinson Micro w ave Anisotropy Pr obe) space 
mission (iNolta et all (|2009| ). iLarson et al (|2010l )'): ACBAR 
(IReichardt et a II I2OO9I ) which has made high resolution mea- 
surements using the Viper telescope at South Pole Sta - 
tion; the South Pole Telescope - SPT (iLueker et alll2009l) : 
the Cosmic Background Imager-CBI jMason et a' 20031 ) 
and the ACT telescope (|Fowier et all I2OI0I ') in the Ata- 
cama desert; the polarisation measurement at Amundsen- 
Scott South Pole Sta tion by the Degree A ngular Scale In- 
terferometer - DASI ijCarlstrom et alll2003l ); as well as bal- 
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Ion bo rn e experiments su ch as BOOM ERANG (jjones et all 
(I2OO6I') .[Masi et all ([200^ ')'). Archeops ("Bcnoit ct all j|2003f ) . 
iTristram et all (120051 ) ) and MAXIMA (Lcc ct al 200l|). Ta- 
ble 1 summarizes typical sky coverage, the accessible £ range, 
and microwave frequency range for some of the above men- 
tioned instruments. 

Determining the CMB angular power spectrum is a 
complex process and the data analysis must take into ac- 
count a number of effects, such as the non-stationnary noise 
contribution, non-circular instrumental beams, the sky scan- 
ning strategy and foreground contamination. A review of the 
general met hods for CMB data proces sing can be found, for 
example, in (|Tristram fc Ganga llioosi '). 

The different systematic and statistical effects on CMB 
angular power specrtum estimates from partial sky maps 
have already been studied in length. In particular, sev- 
eral Montecarlo or analytical methods have been devised 
to evaluate or correct the impact of limite d sky cov- 
erage on estimated angular powe r spec trum (|Hivon et all 
(|2002l ) lMitra et all (|2009l ) lDas et all (|2009l ')). 

In this paper we propose a different approach which is 
based on the analysis of the distortion of the angular correla- 
tion function. However, it should be noted that explicit com- 
putation of the angular correlation function is not needed. 
The use of the correlation function has also been studied by 
several authors, alt hough with a dif f erent a p proach than the 
one dev eloped here (| Szapudi et al I (|2001al 'l , I Szapudi et al I 
(|2001bl )h In section 2, we recall briefly the main relations 
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Table 1 . ^-range and sky coverage for few of the microwave sky observation instruments 



Name 


sky cov. (deg^) 


Angular Resolution (arcmin) 


t range (min-max) 


frequency coverage (GHz) 


WMAP 


~ 40 000 (4 TT sr) 


10-60 


1-1500 


23-94 


ACBAR 


600 (10 fields) 


5 


200-3000 


150,220 


BOOMERANG 


750 


6-10 


10-1500 


145,245,345 


Archeops 


5000 


11 (@143) 


10-700 


143,217,353,545 


CBI 


40 (3 fields) 


5-10 


200-3500 


26-36 


SPT 


100 


~ 1 


2000-9500 


150,220 


ACT 


228 


1.4 


600-8000 


148 



between the angular correlation function ^(7) and the angu- 
lar power spectrum Ct, and express them as a set of linear 
algebraic equations. Using this formalism in section 3, we 
compute the distortion of the estimated Ce in incomplete 
sky maps and we show how the angular correlation function 
can be corrected to minimize these distortions. We show in 
particular that power spectrum estimates on partial maps 
sufTer a loss of resolution in ^-space and we propose near 
optimal ^-space window function (binning). 

The algebraic approach presented in section 3 can not 
be used to compute the variance of the reconstructed angular 
power spectrum (see paragraph 3.1 and 3.4). We have thus 
used Montecarlo simulations to estimate Ce uncertainties. 
We present in section 4 the corresponding results for small 



maps (f2 ~ 47r X 10 ), represe ntative of 



round or balloon 



experiments such as OLIMPO (|Nati et alll2007h . as well as 
the systematic shifts and possible correction for nearly com- 
plete maps (f2 ~ 47r x 0.9) represe ntative of space missions 
such as Planck (|Planck C0U.II2OO6I ). 



2 ANGULAR CORRELATION FUNCTION 
AND POWER SPECTRUM FROM THE 
FULL SPHERE 

We recall here the basic relations for the angular power spec- 
trum and correlation function. Detai led derivation of most 
of th eses formulae can be fo und in (|Magneville fc PansartI 
[2003) or (|WandeIt et aj|200ll ). 

We consider a real signal s(f2) on the sphere (S^), mea- 
sured for each direction Q, — {9,<j)). The function s can be 
expanded on the spherical harmonic basis: 



EE 



aim YP{9, (j>) 



aim — 



(1) 



(2) 



For an isotropic random signal, the angular power spec- 
trum Ct characterizes the statistical properties of the signal 
(() denotes ensemble average). 



{o-lmO-lim') 



(3) 



It is possible to compute an unbiased power spectrum esti- 
mator C( from the spherical harmonic expansion coefHcients. 
The Ci coefficients are independent random variables with 



variance cr~^ (cosmic variance): 



m^+e 



Ce = 



2£ 



— y 

+1 ^ 



{Ci) = Ci 



2 



2 9 
-Cj 



(4) 



(5) 



Ci 2£+l 

For an isotropic signal, the angular correlation function ^(7) 
can also be used to characterize the signal properties, where 
7 is the separation angle (0037 — Q ■ fi') : 



1 °° 

= ^J2m+l)CePi{coa'y) 



(=0 

hi 



(6) 

Ci = 2n J C(7) J'Kcos7) dcos7 (7) 
-1 

It should be noted that the above relation holds also for 
the computed angular correlation function (^(7)) and power 
spectrum (Ce) for a given sky realization. 

00 

= ■^J2(^e + ^)CiPi{cos'y) (8) 
= J s{Q.)s{n')S (cos -f-n.n'^dndQ' 



AA(7) = 



hi 



cos 7 



The above relations, analogous to Fourier series, contain al- 
gebraic (sum) and analytic (integral) expressions, and in- 
volve the discrete variable i" € N, as well as the continuous 
variable 7 £ [0,7r]. We can rewrite these expressions in a 
purely algebraic form, similar to the Discrete Fourier Trans- 
form (DFT) . Indeed, if the Ci spectrum is negligible for large 
£ > £max , and for any set of discrete values {7;} of the 7 
angle, the relations [6] and |8] can be written in matrix form: 



K(7.)] = m = [[Ki.]] * [Ce] 



(9) 



In the particular case where the number of ji values is equal 
to the Imax + 1 non-zero Ce values, the [[Kif]] matrix would 
be square. It can then be shown that the separation angles 7; 
at which the correlation function ^ is computed (^^ — C(7i))i 
could be choosen such that the square matrix [[K]] of size 
i^niax + 1) X {Imax + 1) Is uon-siugular (see appendix). The 
equation may thus be inverted to get: 



[Ci] 



(10) 
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Moreover, these 7^ values are not very different from 
f-max + 1 equidistant 7, distributed from to vr, with the 
matrix elements close to [[K]]^^^ oc Pf(cos7i). 



3 ANGULAR POWER SPECTRUM FROM 
PARTIAL MAPS 

3.1 Angular correlation function distortion 

The above equation (|1UI) can be used to analyse the impact 
of any linear distortion of ^ on the calculated Ci. Any lin- 
ear distortion of the angular correlation function (including 
truncation for 7 > 'ymax) can be represented by a matrix 
[[D]] applied to ^. 

{[^(7)]') = [[D]]*[C(7)1 
{[aY) = [[K]]-' * [[D]] * [[K]] * [Gl 

{[a]") = [[Be]]*[G] 

[[Be]] = [[K]]-^ * [[D]] * [[K]] (11) 

It should be stressed that in most cases, we can only compute 
the mean distortion caused by the measurement process, in 
particular due to the incomplete coverage. The above rela- 
tions will not hold in general for a given realization or a sin- 
gle measurement. It will only be valid when an ensemble av- 
erage is taken for the angular correlation function and power 
spectrum. This explains why it can not be used to estimate 
statistical uncertainties on the computed power spectrum. 

Using the formalism described above, we have computed 
the effect of distorting or modifying the angular correlation 
function in some typical cases such as: 

• Application of a sharp or a smooth cut ^(7) — ^''(7) = 
^(7) *Cut(7) to the angular correlation function. In the case 
of a sharp cut, ^(7) is set to zero for 7 > 7mai (step function 
cut). 

• Undersampling by a factor p > 1, i.e. known only 
for i = k X p 

• Binning effect, where ^(71) is replaced by the mean 
value of ^ in a small interval around 7^. The binning can 
represent the distortion of the ^ when computed through 
the histogram of all pixel pair products s{Q,)s{Q'), binned 
as a function of their separation angle 7 (H • f2' = cos(7)). 

Figure [1] illustrates the effect of restricting the 7 range 
to ^max = 30° {i.e. setting ^(7) = for 7 > 7mai), (left) 
and the undersampling (right). Restricting the range of an- 
gles for ^(7), creates a correlation between different £ while 
undersampling produces an aliasing effect. Both of these ef- 
fects are analogous to well known effects in standard Fourier 
analysis. 



3.2 Partial maps : truncated 5' (7) 

When the angular correlation function is calculated from 
maps with a maximal extent 'ymax, nothing can be known on 
the correlation function for 7 > jmax. In addition, the sta- 
tistical errors on the estimated correlation function will be 
larger compared to the one computed on the corresponding 
full (47r) map. It is well known that computing the angular 
power spectrum from a truncated ^'(7) using the integral 



equation (O or the linear combination equation (|10[) pro- 
duces spurious oscillations. However, these oscillations can 
be filtered out if the resulting power spectrum is binned. 

(dj) = [[Bc]]*[C,] (12) 

Cl = ^«^l(^)C,* (13) 

t 

The filtered or weighted power spectrum Ci defined here 
is obtained by applying the weight function wl(£) to the 
power spectrum Cj. The weight function should be centered 
around L and normalised such that WLii) = 1. The 

wl{(-) would be in general positive, with a maximum for 
£ = L and decreasing to zero ( wl{1) — >■ 0) when | ^ — L | 
increases. The expectation value of the filtered spectrum can 
then be expressed using the weight matrix [[W]]^^ ^ = wl{£): 

{61} = [[W]] * [[Be]] * [C,] (14) 
(Cl) = [[Bcw]]*[Cf] (15) 

We have computed the matrix [[Bcw]] relating true ex- 
pectation values of Ce to mean value of the weighted power 
spectrum {Cj^) from truncated angular auto-correlation 
function. The results shown here correspond to ^(7) trun- 
cated above ^max ~ 30°, for Gaussian or square (step- wise) 
weight functions. Figure [2] shows one of the rows of the 
[[Bcw]] matrix around £ ~ 500 for different weight func- 

TT 

tions. For Gaussian weights width ag > , the effec- 

five ^-space window or filter function becomes numerically 
very close to the corresponding Gaussian function. It can be 
seen also that applying top hat or square weights result in 
an effective window function significantly different from the 
corresponding square function. We have checked that this 
property and the value of optimal weight function width 
ae ~ does not depend on the value of the central L, 

^max 

at least for large enough L > at. 

Our computation and simulations suggest that a Gaus- 

TT 

sian function with a width A£ ^ at — is a near opti- 

mal binning when the angular power spectrum is calculated 
from partial maps with maximum angular extent ~ "imax. 
The word "optimal" should be understood in its common 
sense, and not the mathematical one, as the optimal solu- 
tion depends on the chosen quantitative criteria. For exam- 
ple, different binning should be used if one seeks to increase 
the spectral resolution or if one is concerned by the sta- 
tistical errors. As explained above, the suggested Gaussian 
binning has the following properties: 

• The Gaussian weighted corrected angular spectra re- 
duces the window function tails by an order of magnitude, 
compared to pseudo-Cf , while maintaining the ^-space reso- 
lution as well as similar statistical uncertainties (see section 
|4ll . Using top hat (square) binning yields a window function 
with long tails and oscillations. 

• The resulting true i!-space window function is nearly 
identical to the weight or filter function applied to the com- 
puted Cl, with a simple analytical expression (Gaussian). 
This property is useful for presenting measured power spec- 
tra. 

When the power spectrum is estimated using incom- 
plete spherical maps, a loss of resolution in angular frequen- 
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Figure 1. Color scale representation of the [[Be]] matrix relating {Cj) ^ Cg. Left; Correlation effect induced by restricting the angular 
range to 7 < 30°. Right: aliasing effect when ^ is undcrsampled by factor p=5. The color scale has been chosen to enhance visually the 
matrix structure (the diagonal terms are around ~ 0.15 — 0.2). 



cies (A£ ~ n/'ymax) can also be understood using linear com- 
binations of Legendre polynomial. We define new functions 
Ql, through linear combinations of Legendre polynomials 
around a given L, with Gaussian weights : 

Ql' - I^E,exp[<£^]p. 

- E,exp[ii5^] (16) 

Qj'(cos7) ~ P£(cos7) 7 < 7r/(Tf 

Qp(cos7) ~ 7 > tt/o-j 

Note that these new functions Ql(cos7) behave like 
Pi (cos 7) at small angles (7 < 'ymax ~ t^/o'i) and become 
negligible at large angles (7 > ^max = Tr/ai). This property 
is illustrated on the figures [3] and 2] It might be possible to 
use these Ql functions as a basis to decompose truncated 
angular correlation functions. 

3.3 Partial maps : Decomposition on the Y/" basis 

As we mentioned in section 2, a signal defined as a func- 
tion of direction can be decomposed in spherical harmon- 
ics (see equation {T])). Optimized algorithms for performing 
numerically this decomposition on pixelized spherical maps 
are co m monly used in anal yzing CMB data (|Muciaccia et all 
(|l997h . iGorski et all (|2005h ). The computed aim coefficients 
are used to derive the angular power spectrum Ce using 
equation Q. When an incomplete (or partial) sky map is 
expanded on the Y/^ basis, the resulting power spectrum 
suffers systematic effects and is often called a pseudo— 
spectrum. The pseudo— spectrum distortion has already 
been studied by several authors and can be found for ex- 



ample in (|Hivon et all (|2002l ). lMagneville fc PansartI (|2007h ). 
The partial map can be written as the result of applying a 
mask on the original signal. One can then show that the aim 
coefficients computed on the masked map can be related to 
the true aem through a linear relation, which can be written 
in matrix form: 

s^{e,(l)) = s(6l,<j!.) X mask(e,(^) (17) 
[hem] = [[A]] * [aem] (18) 

The above relation is completely general, independent 
of the statistical properties of the signal s(6,(j>) and valid 
for each realization. The [[A]] matrix coefficients depend on 
the mask spherical harmonic decomposition and Clebsch- 
Gordan coefficients. Then, using the isotropy of the signal, 
it is possible to compute the coupling matrix [[B^^/]] relating 
the pseudo— Cf power spectrum to the true angular power 
spectrum. Although [[A]] is a huge matrix, (~ 2 lO'^ x 2 lO'' 
elements for Imax ~ 1000), only a small fraction of the ele- 
ments contribute to the [[B«/]] matrix. 

[{C!)] = [[B«,l] * [C,,] (19) 

The computations to obtain the [[B^^/]] coupling matrix 
using the above method are rather tedious, while it is possi- 
ble to compute this matrix simply by applying the formalism 
presented in this paper. 

One can construct an estimator ^' of the angular cor- 
relation function using equation ((8| with the integrals lim- 
ited to the partial map wher e the function s(Q.) is mea - 
sured. A lengthy computation (|Magneville fc Pansart|[2007l ) 
shows that this estimator is unbiased for separation angles 
7 < 'ymax- However a fast computation of the angular corre- 
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BcwG matrix row L~500 - Y 

max = 30 d69 



Sum[PI] for ell=20,1 00,300,500 sigma_ell=5 



1 

0.8 
0.6 



WGauss,ellO=502.5 
0_l = 4 
a_l = 6 
O 1 = 8 




460 



480 



500 



520 , 



ell (I) 

Bcw matrix row L~500 - max = 30 d69 



540 



1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 



Bcw row, ell0=502.5 
BcWG Sig_l=6 (circles) 
Gaussian (dastied) 



Bcwsq, Widtli=20.8 (box)] 
Square (dastied) 



V.:-. 



460 



480 



500 



520 540 
ell (I) 



Figure 2. Weighted C^^ calculation truncated ^ at ■^/max = 30°. 
Top: Row of the weighted Bet matrix ([[Bcw]]) around L ~ 500 
for Gaussian weights (-BcWg) with three different widths cr^ = 
4, 6, 8. Bottom: Comparison of BcW matrix row for L ~ 500 with 
the corresponding Gaussian weights (cr^ = 6), in blue and top hat 
(square) weighted (-BcWSgi width = 6 X \/T2 = 20.8) in red. 



PI , Sum[PI] for ellO=100, Slgma_ell=5,10,20 



I ^ 
I 

n 0.001 



0.01 




PI , 1=100 

Sum[PI] 1=100, Sigmaell=5 
Sum[PI] 1=100, Sigmaell=10 

Sum[Pll 1=100. Sigmaell=20 




1e-07 ^ 
1e-08 - 
1e-09 ^ 

1e-10 

0.1 0.2 0.5 1 2 5 10 20 SO 100 200 

Tbeta (deg) 

Figure 3. Lcgondre Polynomial and Qj^ Pi functions, plot- 

ted for L = 100 and = 5, 10, 20. We have plotted |P<?|, |Ql|). 
Notice also that both axes have logarithmic scales. 
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100 


ellO = 
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£ 0.01 
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I 0.001 
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Figure 4. Legendre Polynomial \Pi\ and and the | functions, 
plotted for ae = 5, and £0 = 20, 100, 200, 1500. Notice that both 
axes have logarithmic scales. 

lation function uses equation © or its linear form © and 
the full sphere function should be used. This corresponds to 
equation (|8]) with integrals on the full sphere S'^ x putting 
s(Q) — when Cl is not pointing to the "observed" partial 
map. This estimator (noted ^^), which is of course highly 
biased, will be the one used hereafter. 

As 5' and are related by a multiplication by the cor- 
relation function ^™^'^''(-y) of the mask itself, we can notice 
that the average of the angular correlation function (^'') 
computed on the masked sphere is related to the correla- 
tion function of the complete sphere by; 



i^'i"/)) ~ C(7) (unbiased estimator) 
(^(7)) = 

Here, ^'"^''''(i) can easily be computed using equation ((6)1 
and the optimized algorithms cited below, applied to the 
mask. The distortion matrix in relation can then simply 
be written as a diagonal matrix D"^^^^ which can be used to 
compute the [[B^^/]]. 



D" 



r'*(70 



(20) 



Figure [S] illustrates the distortion of the angular corre- 
lation function estimated from 5^ coefficients computed on 
the masked map. The undistorted angular correlation func- 
tion correponding to a WMAP-like Ci power spectrum, the 
mean value of the distorted correlation function (^) and the 
mask correlation function ^'"^''''(i) are plotted as a function 
of the 7 angle, for a partial 30° x 30° map. 

The analysis of the coupling matrix [[B^^/]] shows that 
the i!-space resolution of pseudo— (7^ spectrum is compatible 
with the optimal resolution A£ ~ tt /'ymax , but has long tails 
or correlation lengths. These long tails are responsible for 
the systematic shifts of pseudo—Ci spectrum relative to the 
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Angular correlation function, partial 30x30 deg maps 



ell-space window function (10=300 , partial map 30x30 deg) 



1 

0.9 
S 0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




Angular Correlation ^(y) 
<^> from CI on masked map 
J;mask(Y) 30x30 deg (dashed line) 




10 



20 



30 40 50 60 

Y - Gamma Angle (degrees) 



Figure 5. Distortion of the angular correlation function when 
computed from spherical harmonic decomposition of a partial 
(masked) spherical map. The partial map has a 30° extension 
on each of the two orthonal directions. Rescaled ^(7) (black) 
and the distorted (^) (red) are shown for a WMAP-like angu- 
lar power spectrum, as well the ratio _ (^)/^™^=''('y) 
(dashed green). Notice that the angular 7 range is limited to 60° 
on the figure. 



true one. This can be seen on figure [6] below, showing a row 
of the coupling matrix Bu' for £ = £0 = 300. 



pseudo-CI 

Cl[ ^cor] = ymax 30 deg 
Cl[ ^cor] = ymax 36 deg 
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ell-space window function (10=300 , partial map 30x30 deg) 
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Cl[ ^cor] = ymax 36 deg 
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350 
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3.4 Partial maps : Correcting pseudo— Cf spectrum 

In some sense, the correlation function ^(7) obtained 
through the pseudo— computation applies a weighting 
function (^™'''''<^) which is inversely proportional to the num- 
ber of available measurement {s{9,(f>)) pairs i.e., to the sta- 
tistical significance of the partial sphere correlation function. 
It is possible to obtain the unbiased correlation function by 
applying the inverse correction ^j^g correlation 

function of the masked sphere, up to a maximum angle 'fmax 
which should be less than the maximum extent of the par- 
tial map. As we increase ^max, the statistical errors on the 
corrected angular correlation function, and the corrected an- 
gular power spectrum increase. There is thus a tradeoff be- 
tween the £-space resolution and the statistical uncertainties 
of the recovered power spectrum Ce- 

The statistical uncertainties (variance) associated with 
the recovered ^ or Ce, depend on the true power spectrum 
and can not be computed using the formalism described 
here. It is however possible to compute the variance for 
specific partial map shapes, such as polar cap maps (see 
iMagneville fc PansartI (|2007l )). In this paper, we have com- 
puted theses uncertainties using Montecarlo simulation and 
the corresponding results are presented in the next section. 
In figure [51 the pseudo— window function is compared to 
the one obtained for the corrected Ce for two values of 'ymax 
(30°, 36°). On partial 30° x 30°, the maximum possible 
value for 'jmax ~ V2 x 30° ~ 42° . Using a correction up to 
7mai = 36° , we obtain statistical errors comparable to those 
associated with the pseudo— Cf, while significantly decreas- 



Figure 6. 30° X 30° partial map : Comparison of B££/ matrix 
row for £ = 300 (black) with Gaussian weighted C' {£) estimated 
from §*(7), truncated at ■ymax = 30° (magenta) and fmax = 36° 
(Top: linear Y-scale, bottom: logarithmic Y-scalc). 



ing the systematic effects on the recovered power spectrum 
(See section [J] ). 

We propose the following procedure to recover an unbi- 
ased or corrected power spectrum, using fast spherical har- 
monic decomposition of masked spherical maps. 

(i) Compute the pseudo— angular power spectrum on 
the masked sphere through fast spherical harmonic de- 
composition, as well as the pseudo— (7^ of the the mask itself. 

(ii) Compute the corresponding discrete angular correla- 
tion function and the mask correlation function 



mask 



'Cf 



(21) 



(iii) Define the truncation angle -ymax compatible with the 
maximum map extent and desired £ resolution. Compute the 
corrected-truncated angular correlation function using the 
diagonal correction matrix Door. 



(22) 



[e(70] = [[Dcor]] * [e''(7, 
Dcor(i,i) = l/r"*(70 7. 
Dcor(i,i) = i ^ j or "n > ^max 

The correction matrix is equal to the inverse of the mask 
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distortion (see equation (|20|l ') if angular correlation infor- 
mation is present up to jmax ~ tt. It is useful to apply a 
step smoothing function to avoid the discontinuity of ^' for 
7 = Jmax, such as a sigmoid: 



/(7) = 1/ [ 



1 + exp 



7-7" 



with S ~ O.OI7 rr 



This step smoothing function enhances the behavior of the 
computed power spectrum, decreasing residual oscillations. 

(iv) Compute the corrected angular spectrum and apply 
the Gaussian filter function in ^-space with at = n/ymax- 



K 



Wg(4,^) 



[C'(70 



(23) 



[[Wgi: 

^ exp((<!-<'o)V2,T|) 
Ef cxp((^-<'o)2/2ct|) 

The coupling matrix Be is independent of the true power 
spectrum and can be computed using the following relation: 



{[Cl\) = [[Bc]]*[C,] 
[[Be]] = [[Wg]]*[[K- 



IDcor]] * [[K]] 



(24) 



4 SIMULATION RESULTS ON SMALL MAPS 
AND MASKED MAPS 

There is a tradeoff between the achievable resolution and 
the uncertainties of the estimated power spectrum. In order 
to evaluate these uncertainties, we have performed Monte- 
carlo simulations to generate partial and full sky spherical 
maps and compute the power spectrum using spherical har- 
monic decomposition and the corrected/truncated angular 
auto correlation function. The method to correct system- 
atic effects proposed in this work is independent of the true 
power spectrum. The computation of the coupling matrix 
and correction matrix is rather fast with a CPU time com- 
parable to few Monte-Carlo realisations using fast spherical 
harmonics map analysis (~ minute). Reasonable estimates 
of errors for a given input power spectrum can be obtained 
using few hundred realisations. Although it is possible to 
compute the systematic effects due to the limited sky cover- 
age using Monte-Carlo, this would require much larger com- 
putation time compared to the uncertainty estimates. The 
estimate of the [[B^^/]] coupling matrix elements with suffi- 
cient precision would need a large number of Monte-Carlo 
realisations, due to the cosmic variance. In addition, this 
computation process would have to be repeated for several 
input power spectra. 

The simulations have been performed with different par- 
tial map geometries and input (true) Ce power spectra to 
check the validity of the conclusions. However, for the sake 
of clarity, we present here the results only for two map ge- 
ometries and two input spectra. We have used WMAP-like 
true Ce power spectrum, labeled tmmap and a second shape, 
labeled te, which have sharp features on top of a smooth 
spectrum C{£) = Cte/[£ x (^ -|- 1)] -I- peaks in order to illus- 
trate the ^-space resolution effect. The results presented here 
have been obtained for the following two map geometries: 

(i) Partial square map, 30° x 30°, covering the angular 
range 75° ^ 9 ^ 105° and 00 sS (ji 0o + 30° . This map cov- 
ers ~ 0.27 sr — 0.022 x 4n sr. This small map corresponds 
to the case of ground or balloon CMB instruments. 
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Figure 7. WMAP KQ85 CMB temperature analysis mask (top) 
and the corresponding normalised correlation function gmask^.^^ 
(bottom). 



(ii) A nearly complete spherical map, but with an equa- 
torial band removed (set to zero). This case shows the dis- 
tortion of the C{£) power spectrum obtained from all sky 
space experiments such as WMAP or Planck, where some 
part of the sky with strong non-CMB microwave signals has 
to be masked. The removed equatorial band illustrates the 
effect of Galactic cut applied to CMB maps. We have per- 
formed simulations on 47r spherical maps with a — ±10° 
equatorial band cut (s(6l,0) = for 80° < 6* < 100° ,V<;/)). 
Such a map covers 82.5 % of the whole sky (0.825 x 47r sr) 
and provides correlation information for the whole angular 
range (0 < 7 < vr). 

(iii) A more realistic case where we have applied the 
5-year WMAP C MB temperature analysis (KQ85) mask 
(|Gold et alll2008h to simulated CMB maps. The mask and 
corresponding normalised correlation function is shown in 
figure [7] The unmasked area corresponds to ~ 0.82 x 47r sr, 
nearly equal to the area of maps with a ±10° equatorial 
band removed. It has a complex and patchy shape. 

Given the overall shape of the Ci spectra studied here, 
we have represented £x (^+1) Ce for the power spectrum, and 
£ X (£+ 1) ace for the associated statistical uncertainties on 
all figures that follow in this section. The systematic shifts 
of the recovered power spectrum can be quantified using 
difference between the recovered and true power spectrum, 
normalised to the cosmic variance acv{£) 
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Figure 8. 30° X 30° partial maps : Comparison of computed 
power spectrum and true t^map power spectrum (top) and the 
associated statistical errors (bottom). True power spectrum in 
black, pseudo— in red, C| computed from corrected angular 
correlation function truncated at 'imax = 30° (light blue), and 
corrected C{lY binned with Gaussian weights (blue) 

4.1 Small 30° x 30° maps 

The figures [S] and [5] show the recovered power spectrum on 
small 30° x 30° maps for a true WMAP-like power spectrum. 
The systematic shifts of the pseudo— power spectrum are 
clearly visible, in particular on figure[5](top). The false oscil- 
lations on the raw power spectrum recovered from truncated 
^ is also shown on the figure [5] (top). 

It can also be seen that by choosing "fmax = 36°, it 
is possible to avoid nearly all the systematic shifts of the 
power spectrum, while keeping the level of statistical fluc- 
tuations comparable to the pseudo— Cf , and without losing 
the ^-space resolution (See figure |9] and [10}. However, it 
should also be noted that the corrected Ci suffers a system- 
atic overestimation at low I < (2 — 3xo-f = 2 — "iir / 'jmax) , 
as expected from the coupling matrix Be. 

The systematic shift 6{t) decreases by a factor 10-30, 
changing from 5 ~ 0.25 — 0.75 in the case of pseudo— to 
less than 0.05 for Gaussian weighted C\ computed from ^(7) 
corrected up to 'ymax = 36° (figure [9]). 

4.2 Maps with ± 10° Galactic cut 

The figure [TT] shows the bias introduced on the recovered 
power spectrum, specially at low multipoles ^ < 30 using un- 
corrected C\ from map projection on the basis. In this 
case, the angular correlation function can be estimated at all 
angular scales, but the ^(7) obtained from the pseudo— is 
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Figure 9. 30° X 30° partial maps : Comparison of computed 
power spectrum and true twmap power spectrum (top) and the 
associated statistical errors (bottom). True power spectrum in 
black, pseudo— in red, Gaussian weighted binned Cf computed 
from corrected angular correlation function truncated at "fmax = 
30° (blue), and with ■ymax = 36° (violet) 

distorted. By correcting the ^(7), it is possible to recover the 
unbiased power spectrum Cf, as shown on figure [TT1 (top). 
The systematic shift S{£) changes from S ~ 0.25 in the case of 
pseudo— to less than 0.05 for Gaussian weighted Ce com- 
puted from corrected ^(7). However, the statistical fluctua- 
tions are larger, compared to the cosmic variance, reachable 
if the complete An map is available, or to the pseudo— 
variance on the cut map. Figure [121 shows the effect of the 
equatorial cut on a sharp feature present in the simulated 
power spectrum around £ = 50. The power spectrum com- 
puted from corrected angular autocorrelation also corrects 
this distortion. 

Although more sophisticated methods can be used to 
recover the spectrum at low multipoles from cut maps, the 
correction algorithm proposed here can be used to easily 
correct the pseudo— Cf spectrum. 

4.3 Maps with WMAP KQ85 mask 

The effect of the WMAP KQ85 foreground suppression mask 
on the recovered power spectrum is shown in flgure [T3l This 
mask cuts an equatorial area on the map which is signifi- 
cantly less extended than in the case of the simple ± 10° 
cut discussed above, but it removes also a large number of 
smaller patches of the sky scattered over the map. As a 
result, it can be seen that the systematic effects at low mul- 
tipoles on the recovered pseudo— are smaller compared 
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Figure 10. 30° X 30° partial maps : Comparison of computed 
power spectrum and true tg power spectrum (top) and the asso- 
ciated statistical errors (bottom). True power spectrum in black, 
pseudo— Cf in red, Gaussian weighted binned C| computed from 
corrected angular correlation function truncated at "fmax = 30° 
(blue), and with 'ymax = 36° (violet) 

to the more extended ± 10° equatorial cut, while significant 
shift appears at higher multipoles, around £ = 400 for in- 
stance, due to the patchy structure of the mask. This patchy 
structure is also responsible for the oscillatory behaviour of 
the Cf power spectrum computed from the corrected angu- 
lar correlation function, for i > 10. Theses unphysical oscil- 
lations can be smoothed out by applying a narrow Gaussian 
{ere = 0.75) filter function, C^'^'^ , as can be seen on figure [T3l 

As expected, correcting ^(7) using ^™^'^'^('y) increases 
the statistical uncertainties compared to pseudo— C<, while 
the filtered power spectrum has smaller variance. It is pos- 
sible to use the combined power spectrum, C|'^ without fil- 
tering for < 10 and Gaussian weighted C^'^''^ for £ > 10. 
The systematic shift S(i) decreases from 5~0.5— 1.5in the 
case of pseudo— to (5 < 0.05 for this combined Ce- 

We have also checked that our conclusions are valid in 
the presence of noise. We have performed simulations where 
Gaussian fiuctuations, with smooth spatial variations have 
been added to simulated maps. The reconstructed power 
spectrum corresponds to the sum of the sky signal and the 
noise spectra {Ce = Ce + C"™"''), distorted by the mask, as 
it is expected for nearly isotropic uncorrelated noise. 

5 CONCLUSIONS 

We have established a simple method to evaluate and correct 
systematic effects associated with power spectrum computed 
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Figure 11. Att maps with it 10° equatorial cut: Comparison 
of computed power spectrum and true t^map power spectrum 
(top) and the associated statistical errors (bottom). True power 
spectrum in black, pseudo— in red, Gaussian weighted binned 
C^"^ computed from corrected angular correlation function in blue. 



using the pseudo— C7f on partial maps. The coupling matrix 
Bee' can indeed be computed using matrix algebra and the 
sky mask spherical harmonic decomposition. We have also 
described an algorithm to correct the systematic shifts of 
the calculated power spectrum, as well as a near optimal 
^-space window or filter function. It should be noted that 
it is possible to improve the I resolution, compared to the 
natural resolution /S.I ^ ai ^ it /'ymax at the expense of 
higher statistical fluctuations. For all sky CMB experiments 
such as WMAP or Planck, non-CMB dominated parts of 
the sky (galaxy . . . ) are usually excluded from the angular 
power spectrum estimation. We show that our method can 
be used to correct the power spectrum distortions at low I 
in such cases, when the angular power is computed using 
fast spherical harmonic decomposition on almost complete 
maps. 

The correction method described here can easily be ex- 
tended to the CMB-polarisation power spectrum. A sim- 
ilar approach for polarised maps has been developed in 
(|Chon et all 1200 j ). It should also be possible to use this 
method to improve the angular power spectrum calculation 
by taking into account individual pixel measurement errors. 
For observations with negligible correlated noise, a weight- 
ing function mask inversely proportional to the individual 
pixel measurement uncertainties can be applied to the map 
before decomposition on the Yem basis, and subsequently 
corrected for. 
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Figure 12. Air maps with it 10° equatorial cut: Comparison of 
computed power spectrum and true tg (C(£) = Cte/l X + 1) + 
peaks) power spectrum (top) and the associated statistical errors 
(bottom). True power spectrum in black, pseudo— in red, C^'^ 
computed from corrected angular correlation function in blue 



APPENDIX 

We will show that one can invert equation Q for a certain 
choice of the separation angles 7. We assume that the power 
at large £ is negligible. We can then consider that Ci = Q for 
all multipoles with I > Imax- The sum in equation ((6]) be- 
comes finite, so for given discrete values of 7, the matrix 
elements in equation @ are: 



with 1 ^ i and ^ ^ ^ imax- 

If the matrix [[K]] is invertible, one may write [Ci] — 
[[K]]~^ * l^i] ^-^id the Ct spectrum can be recovered from 
the values of the angular correlation function computed at 
angular separations 7^. 

In the following we will show that the separation angles 7^ 
can be choosen such that the square matrix is non-singular. 

Pi{x = cos (7)) is a polynomial of degree £ in the interval 
[—1, Using equation ^ we see that ^(7) is a polynomial 
of degree at most imax ■ The integrand ^ x of equation ((?]) 
is a polynomial of degree at most imax + 1, so for the highest 
multipole Ci^^^ the integrand has a degree at most 2£max- 

We know that, using Gauss-Legendre quadrature of or- 
der n, integrals of polynomials of degree up to 2n -j- 1 can be 
exactly expressed as a n-term weighted s um of the polyno- 
mials computed at special values of x (see lAbramowitz et aj 

|i9z3)). 

For all £ ^ £max and for L ^ £max + 1, the integral in 
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Figure 13. 47r maps masked with the WMAP 5 year tempera- 
ture KQ85 mask. Comparison of recovered power spectrum and 
true tuijnap power spectrum (top and middle) and the associ- 
ated statistical errors (bottom). True power spectrum in black, 
pscudo— in red, dashed, C|'^ computed from corrected angular 
correlation function in blue, and filtered with Gaussian weights 
{ae = 0.75) in violet C^"'^ . 



equation (O can be expressed as the algebraic sum: 



(26) 



where Xi = cos(7i) are the L roots of Pl{x) and Wi the 
weights. For L ^ £m.ax + 1, this can be written in a matrix 
form: 



K'li = 27rPf (cos(7i)) x Wi 



^ ^ ^ ^max 



Clearly, if we take L = £max + 1, the number of abscissas 
values is n-y — £max + 1 and the matrix \[K]\ is square. In 
that case [[A'']] = [[A']]^^ and equation ((Tjl can be written 

[Ct] = [[A]]-' * [e,] (1 < i < £max -t 1, < £ < £max) 

(27) 

Remark: Numerically, the angular correlation function 
could also be computed at £max + 1 regularly spaced val- 
ues 7i between and tt. This is possible because the 7; cor- 
responding to the Legendre polynomial roots Xi = cos (7;) 
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satisfy the relation This paper has been typeset from a T^^X/ ffl^^X file pr 

2i - f T 2i + f by the author. 

and are nearly regularly spaced. 

Thus the angular distance between a Legendre polyno- 
mial root and the nearest regularly spaced abscissa is always 
lower than the resolution associated with the highest multi- 
pole (~ -K/tmax)- 
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